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ABSTRACT 

Accumulating evidence indicates that some of ultra-luminous X-ray sources (ULXs) are intermediate mass 
black holes (IMBHs), but the formation process of IMBHs is unknown. One possibility is that they were formed 
as remnants of population III (Pop III) stars, but it has been thought that the probability of being an ULX is too 
low for IMBHs distributed in galactic haloes to account for the observed number of ULXs. Here we show that 
the number of ULXs can be explained by such halo IMBHs passing through a dense molecular cloud, if Pop 
III star formation is very efficient as recently suggested by the excess of the cosmic near-infrared background 
radiation that cannot be accounted for by normal galaxy populations. We calculate the luminosity function 
of X-ray sources in our scenario and find that it is consistent with observed data. Our scenario can explain 
that ULXs are preferentially found at outskirts of large gas concentrations in star forming regions. A few 
important physical effects are pointed out and discussed, including gas dynamical friction, radiative efficiency 
of accretion flow, and radiative feedback to ambient medium. ULXs could last for ^ 10^~* yr to emit a total 
energy of ^ lO^'' erg, which is sufficient to power the ionized expanding nebulae found by optical observations. 
Subject headings: black hole physics — X-rays: galaxies — infrared:general 



L INTRODUCTION 

Ultra-Luminous X-ray sources (ULXs, Makishima et al. 
2000 and references therein) are bright X-ray sources hav- 
ing luminosities greater than ^ 3 x lO-'^ ergs s"' found in off- 
nuclear regions of nearby galaxies. The luminosity exceeds 
the Eddington limit of a 20Mq black hole (BH) that is 
the maximum mass expected from normal stellar evolutionary 
paths (Fryer & Kalogera 2001), and their origin is now a mat- 
ter of hot debate. A few models have been proposed to explain 
the high luminosity from ordinary stellar mass black holes, 
such as geometrical beaming (King et al. 2001; King 2004) or 
super-Eddington accretion (Begelman 2002). However, clear 
observational evidence for these has not yet been found. On 
the other hand, the high luminosity may also be explained by 
intermediate mass black holes (IMBHs) of Mbh ^ IOOM0, 
without violating the Eddington limit. Recently evidence for 
IMBHs has been accumulating for many ULXs. They show 
lower temperature (^ 0. 1 keV) blackbody components than 
that of the Galactic black hole binaries, which is consistent 
with IMBHs being ULXs (Miller et al. 2004). It is known 
that, for X-ray binaries and active galactic nuclei (AGNs), the 
break frequency of the power density spectrum of X-ray flux 
variability scales with the black hole mass, and the break fre- 
quency of ULXs is indeed between those of black hole bi- 
naries and AGNs, again indicating intermediate BH mass for 
ULXs (Cropper et al. 2004). It appears that at least a part of 
ULXs, especially the most luminous ones, are IMBHs. (See 
e.g., Colbert & Miller 2004; Mushotzky 2004 for recent re- 
views.) 

However, the formation of IMBHs to become ULXs is a 
challenging problem. One possibility is formation in dense 
stellar regions or stellar clusters (Portegies Zwart & McMillan 
2002). However, though ULXs are preferentially found in 
actively star forming regions, their position is near but not 
coincident with such stellar clusters, disfavoring this scenario 

Electronic address: mii@kusas tro.kyoto-u.ac.jp1 



(Zezas & Fabbiano 2002). 

On the other hand, it is theoretically expected that IMBHs 
form as remnants of massive population III (Pop III) stars 
(e.g., Schneider et al. 2002), and they may become ULXs if 
they accrete sufficient gas in nearby galaxies (Madau & Rees 
2001). However, it is not easy to explain the observed number 
of ULXs by this scenario. Islam, Taylor & Silk (2004) con- 
sidered a possibility that free-floating IMBHs accrete from a 
residual core of gas associated with the IMBHs from their for- 
mation epoch. It is highly uncertain whether such gas core can 
survive until present from the early universe, and this scenario 
cannot explain the majority of ULXs that are associated with 
star forming regions. Volonteri & Perna (2005) recently per- 
formed a similar calculation, and they found even more pes- 
simistic prediction for the observable ULX number Krolik 
(2004) considered a possibility that free-floating IMBHs be- 
come ULXs when they enter into dense molecular clouds in 
star forming regions of galactic disks, so that the association 
with star forming regions can be explained. Krolik consid- 
ered a hypothetical "disk" population of Pop III stars that form 
within a galactic disk, from primordial gas that has somehow 
been preserved until the disk formation. This is because the 
number of ULXs would be too small if, as normally expected. 
Pop III remnants are distributed in a galactic halo. It is an 
interesting possibility, but it is difficult to quantitatively esti- 
mate how much amount of such gas can survive during the 
formation of galactic disks. 

The efficiency of Pop III star formation is highly uncertain, 
but there are some observational hints. Independent groups 
(Wright & Reese 2000; Cambresy et al. 2001; Matsumoto 
et al. 2004) reported detections of the cosmic near-infrared 
background radiation (CNIB), which cannot be explained by 
normal galaxy populations (Totani et al. 2001). The excess 
is particularly large at ^ 1.5/im, which is at ~ 5-6(7 level ac- 
cording to the IRTS data of Matsumoto et al., though there 
may still be large systematic uncertainty in the CNIB mea- 
surements, such as zodiacal light subtraction. It has been 
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shown that efficient formation of Pop III stars before z ^ 8 can 
explan this unaccounted excess of CNIB (Santos, Bromm, & 
Kamionkowski 2002; Salvaterra & Ferrara 2003); the spectral 
peak of the excess at ^ l.Sfim corresponds to the redshifted 
Lyman a line emission. 

The required amount of Pop III stars to explain the CNIB 
is somewhat extreme (see, e.g., Madau & Silk 2005); about 
~ 10% of the cosmic baryons must be converted into Pop III 
stars. However, paucity of theoretical understainding of Pop 
III stars does not allow to reject this possiblity completely, 
and it is the only one realistic explanation of the CNIB ex- 
cess, except for the systematic uncertainty in the CNIB mea- 
surements. Theoretical studies on Pop III star formation in- 
dicate that they are likely more massive than normal stars, 
and most of Pop III stars could be massive enough to collapse 
into IMBHs without any metal ejection, meeting the metal- 
licity constraints from quasar absorption line systems. The 
consequence of this scenario is that there should be IMBHs 
with Mbh ~ 1O^-1O"^M0 in halos of nearby galaxies, whose 
mass fraction in the baryonic matter is /m = r^iMBH /^b ^O.l, 
where fls ^ 0.05 is the cosmic baryon density in the stan- 
dard definition. This is consistent with observational con- 
straints from gravitational lensing experiments (Wambsganss 
2002) and from stellar dynamics such as heating of disk stars 
(Wasserman & Salpeter 1994), disruption of globular clusters 
(Murali et al. 2000), and disruption of halo wide binaries (Yoo 
et al. 2004). 

In this paper we show that, if this is the case, such IMBHs 
passing through dense molecular clouds of star forming re- 
gions should have a siginificant contribution to the observed 
ULX population, especially to the highest luminosity ones. In 
^ we give an order-of-magnitude estimate of the expected 
ULX number. We discuss two relevant physical processes in 
^(gas dynamical friction) and ^^(radiative efficiency of ac- 
cretion disks). We will present a more detailed calculation 
of the predicted ULX luminosity function in ^ Radiative 
feedback and the ULX duty cycle will be discussed in ^ Fi- 
nally some discussions, conclusions, and predictions will be 
presented in ^ 

2. NUMBER OF ULXs: AN ORDER-OF-MAGNITUDE ANALYSIS 

Detectability of isolated BHs accreting interstellar gas has 
been studied in a number of publications (Fujita et al. 1998, 
Agol & Kamionkowski 2002 and references therein), and we 
follow these for the basic formulations. As a typical galaxy, 
we choose parameters of our Galaxy with a total halo mass 
of Mhaio = 1O'^M0 and circular velocity of Vc = 220 km s~'. 
Luminosity of a black hole with mass Mbh expected from the 
Bondi-Hoyle accretion in a gas cloud is given by 

Lx(n,v) = mc rj= — 5 

= 4 X 10^^r;_i Ml 112 ergs s"' , (1) 
where m is the gas mass accretion rate, v the IMBH veloc- 
ity relative to the gas cloud, and p, auc, and Cj are the den- 
sity, the turbulent velocity, and the sound speed of a molecu- 
lar cloud, respectively. In the third equation, we used M2 = 
Mbh/(100Mo), no = n/(100cm-3), and vx = €;/(10km s"'), 
where n = p/{p.mp) is the hydrogen number density, lUp the 
proton mass, /it = 0.13 the mean molecular weight, and v = 
(i;^ + cr^(^ + cj)'/^. The parameter 77 = O.I77-1 represents not 
only the radiative efficiency of accreting gas but also in- 
cludes the uncertainty in the Bondi-Hoyle formula of accre- 
tion rate. The turbulent velocity ctmc scales with density as 



(Tmc — 'i-Tn2 '^^ km s ' (Larson 1981). The sound speed c, 
depends on the gas temperature, which is Cj ^ 0.3 km s"' for 
cold molecular gas {T = 10 K) but ^10 km s"' for ionized hot 
gas {T - 10'* K). 

The number of ULXs in the Galaxy, whose luminosity and 
ambient gas density are Lx and n, respectively, is then given 



a'jVuLx 

dndLx 



IMBH /disk 



df„{n) df^jv) 
dn dv 



dv(n,Lx) 



dLx 



, (2) 



where Mmbh = fni(^B/^M)Mha\o/MBH is the total number of 
IMBHs in our Galactic halo, the cosmological density pa- 
rameter Hm = 0.3, /disk the fraction of halo IMBHs passing 
through the disk region, df„/dn the volume fraction of re- 
gions in the Galactic disk which is occupied by molecular 
clouds of density n per unit cloud density, and df^/dv the 
velocity distribution function of IMBHs. In this equation, v is 
a function of n and Lx as determined by eq. Q. 

Assuming the Navarro, Frenk, & White (1997) density pro- 
file with parameters given in Klypin, Zhao, & Somerville 
(2002), we estimate the dark matter mass included within 10 
kpc radius is about 4% of Mhaio ' ■ Then /disk can be estimated 
as - 0.04 X (3/4) 2Hg/(lO kpc) - 4.5 x lO"'*, where Hg = 
75 pc (Sanders, Solomon, & Scoville 1984) is the disk scale 
height of molecular gas. The velocity of ULXs must be much 
smaller than the dispersion of the Maxwell distribution, asH = 
Vc/V2 = 160km/s, and hence t//^,/t/i; - (2/7r)i/V/cr^^. 

Following Agol & Kamionkowski (2002), we model /, as 
dfn/dn cx rT^ with /3 = 2.8, between Mmin =10^ and Wmax = 
10^ cm"-'. Normalization is given as: 
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where (Smc) is the mean surface mass density of molecular 
gas in the Galactic disk, which is averaged for IMBHs passing 
through the disk as: 

/ Iwr Y^Mcir) Pdm dr 



J 2nr Pdm dr 



(4) 



where r is the galactrocentric radius and p^M the dark halo 
density. By using the data of Sanders et al. (1984) for Emc('"); 
we found (Smc) = 29M0 pc"^. 

We cannot take all value of n between and Wmax, be- 
cause V is limited as {) > (o'mc"'"''?)'''^ ^'^'^ " must be larger 
than 



; 2.5 X WL3C)T]Z\Mf 



(lOkms-i) 

where 



(5) 



for v(n,Lx) to have a solution 
Lx/{IQ^'^ ergs s"'). For n greater than n'^^ 



v(n,Lx) - 3.4 ri/^Mll^n^^L^^'^ km s"* 



-39 



and hence we find 



dn dLx 



cx n'-^Lf 



(6) 



(7) 



Therefore integrating over n between Mmin and Wmax, where 
nmin = max(ninin, «[„;„), we find the expected total number of 

' The infall of baryonic gas into the center of a dark halo leads to adiabatic 
compression and a subsequent enhancement of dark matter density including 
halo IMBHs, which will be discussed in ^ 
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ULX as: 

dN 

NvLxO Lx) ~ —Lx (8) 
dLx 

„ Wmin dfn(nmm) /• r , t \i /hn 

^MMBH/disk-^-^ J fv [< U(nmin,ix)] (9) 

p — 2 dn 

= l.lxlO-2L3lry_lM2f^)'^ (10) 

\ ^min / 

If ((t|^(- + Cj)'''^ ~ 10 km s"', n^;^ is larger than nmin and the 
number is reduced to A^j/zjc ^ 8.6 x 10""* fortheMgH = IOOM0 
case, but n^;„ < and no reduction of A^ulx for Mbh = 
IO^Mq. 

There is no ULX in our Galaxy, and the number is known 
to be proportional to star formation rate (SFR) of galaxies 
(Grimm, Gilfanov, & Sunyaev 2003). Under a reasonable as- 
sumption that the amount of molecular gas clouds in a galaxy 
is proportional to SFR of the galaxy, the number of ULX 
should also scale with SFR^. Then /„ of a starburst galaxy 
like the Antennae is 30 times larger than that of our Galaxy. 
The number of ULXs in such a galaxy is then estimated to be 
--3.3 for Mbh = lO^M©, while in the Antennae ^ 10 ULXs 
are observed above lO-''^ ergs s~' (Zezas & Fabbiano 2002). 
This is an encouraging result; it should be noted that the ob- 
served ULXs probably include ordinary X-ray binaries in the 
luminosity range of ~(l-3)x 10'''' ergs s~'. In the following 
sections, we discuss two more physical processes; one (gas 
dynamical friction) may enhance the ULX number and the 
other (radiative efficiency of accretion) is important to repro- 
duce the observed luminosity function of X-ray sources. 

3. THE EFFECT OF GAS DYNAMICAL FRICTION 

The expected number of ULXs is sensitively dependent on 
the relative velocity v, since Lx oc and fv{< v) cx v^. 
Therefore, if there is a process that reduces the speed of 
IMBHs in molecular clouds, it would significantly enhance 
the expected ULX number. Dynamical friction is a possible 
effect, and considering that the star formation efficiency in 
molecular clouds is generally low (stellar to gas mass ratio 
Mstai/MMc ^ 0.01-0.1, Myers et al. 1986; Larson 1988), the 
dynamical friction by gas will be more efficient than by stars. 
The gas dynamical friction time scale is given by (Ostriker 
1999) 

fj. - " - ^^^^ - (11) 

V Fdf 4Tr^G'^MBHpln(rtn-ax/rmin) 

= 8.5 X 10**r't^i^2'n2' yr , (12) 
where F^f is the friction force, vi = ^/(lOkm s"'), r^ax and 
r^iin are the sizes of surrounding medium and the object re- 
ceiving the force, respectively, and ^ is a parameter of order 
unity to take account of the uncertainty in the formula. Here 
we took Tmax = 20 pc and rmin =10^ km Schwartzschild 
radius-') for the Coulomb logarithm, ln(rmax/'"miii) ~ 27. 

This should be compared to the crossing time of molecular 
clouds, 

fci-oss ^ ^ = 2.4 X lOV''"!' yr , (13) 

V 

^ When the morphology of a starburst galaxy is greatly different from that 
of our Galaxy, a simple scaling of the ULX number with SFR may not be 
exactly correct. However, it is difficult to quantitatively discuss this effect 
and it is beyond the scope of the paper. 

Gas within the accretion radius r^i ~ GMbh/v- ~ 10'* cm is captured 
by the IMBH, and momentum loss by the accretion leads to a similar formula 
that is different only by numerical factors. 



where the size of typical molecular clouds is given as Lmc ~ 
24.7«2" '' pc (Larson 1981). The ratio 

-^ = 350r't^f«2°'A^2' (14) 

Across 

becomes of order unity for Mbh ~ 1O^M0 and v <5 km s"' . 
The dependence on v and n indicates that this is achieved 
only for objects having the highest accretion rate, i.e., ULXs. 
Therefore this process could boost up the ULX luminosity 
particularly at the brightest end of the ULX luminosity func- 
tion. 

4. RADIATIVE EFFICIENCY AND LUMINOSITY FUNCTION SLOPE 

So far we assumed that the X-ray luminosity is propor- 
tional to the Bondi accretion rate, i.e., a constant radiative 
efficiency, rj. However, in this case we have a serious prob- 
lem about the ULX luminosity function. Since oc m cx v~^, 
while /u(< w) cx u^, a robust prediction is that the luminosity 
function should scale as dN/dL cx L"" with a = 2. On the 
other hand, the observed luminosity function of ULX in star 
forming galaxies is smoothly connected"* to that of high mass 
X-ray binaries (HMXBs), with a slope of a = 1 .6 in a wide 
luminosity range of Lx lO^^-lO^** ergs s"' (Grimm et al. 
2003). Therefore, even if we successfully explained the ob- 
served number of ULXs, it would overproduce the observable 
number of low-luminosity X-ray sources by the same IMBH 
population. 

This problem can be solved if we take into account the 
change of accretion mode at low accretion rate compared with 
the Eddington rate; recent studies established that accretion 
disks change their state from standard thin disks to radiatively 
inefficient accretion flows (RIAFs), at luminosity below ^ 
10% of the Eddington luminosity (Kato, Fukue, & Mineshige 
1998; Narayan 2004). In the RIAF regime, the luminosity 
scales roughly L cx m^, and hence the luminosity function 
slope of ULXs becomes a = 1 .5, avoiding the overproduction 
problem of low luminosity sources. 

5. MODELING THE LUMINOSITY FUNCTION 

Here we calculate the ULX luminosity function by integrat- 
ing numerically over «, v, and the location of a BH in the 
Galaxy. The calculation is based on the formulation of Agol 
& Kamionkowski (2002), but we also include the follwing 
two effects discussed above, (i) To make the luminosity pro- 
portional to m~ in the RIAF regime, the radiative efficiency 
is reduced as 7] = m/(0.1mEdd)'7Edd when m < O.lmEdd, while 
V = rjEdd for m > O.lmEdd, where mEdd = iEdd/(c^?7Edd), ^Edd 
the Eddington luminosity, and 77Edd is a constant parameter 
to determine the efficiency around the Eddington accretion 
rate, (ii) The relation v(n,Lx) is replaced by v'{n,Lx) to take 
into account the dynamical friction, where v' and v (v' > v) 
are the velocity before and after an IMBH passes through a 
molecular cloud having the size Lmc (see whose rela- 
tion is determined by solving the equation of motion with the 
dynamical friction. Since the initial mass function (IMF) of 
Pop 111 stars is highly uncertain, we simply calculate for a 
(5-function IMF with Mbh = 10^ and lO^M©. 

This fact is often used to argue that ULXs are physically an extension of 
normal HMXBs (e.g., Gilfanov 2004), but the latest work on ULX luminos- 
ity function indicates that ULXs with Lx ^ lO"*" ergs s"' might be a different 
population from lower luminosity ones (Liu, Bregman, & Irwin 2005). In our 
scenario, the similarity between the number of IMBH-ULXs and extralopa- 
tion of the HMXB luminosity function must be a coincidence. 
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The results are shown in Fig.^ We scale up the luminos- 
ity function of the Galaxy to that for a galaxy with SFR>5 = 
5OM0 yr"', to compare with the observed data of the univer- 
sal luminosity function constructed by X-ray sources in dif- 
ferent galaxies and normalized by SFR, where SFR>5 is SFR 
for stars heavier than SMq and SFR>5 = 0.25Mq yr"' for our 
Galaxy (Grimm et al. 2003). Considering the uncertainties, 
we show eight calculations with ^ = 1 or 10, c, = 0.3 km s~' 
or lOkms"', and rj^^^ = 0.1 or 1. It should be noted that 
the parameter rjEdd also includes the uncertainty in the Bondi- 
Hoyle accretion formula. Still, the case of r/Edd = 1 is rather 
unlikely, but we take this value as a maximally possible one. 
The choice of ^ = 10 may also be extreme, but we took this 
value to show the effect of dynamical friction clearly. There 
are two breaks in the model curves: one corresponding to the 
transition luminosity between the two accretion modes (the 
standard thin disk or RIAF) at Lx = O.lLsdd, and another cor- 
responding to the luminosity above which n'^^^ > n^in- As ex- 
pected, the dynamical friction effect boosts up the luminosity 
function at the brightest luminosity range. 

These results show that at least a part of the brightest 
ULXs (Lx > lO'*" ergs s"') can be explained by our scenario, 
without overproducing low luminosity sources in the range 
of Lx ^ 10^^ ergs s"\ where normal HMXBs are dominant. 
(Therefore we do not have to explain all of the low luminos- 
ity sources.) Especially, with Mbh = IO^Mq, the number of 
the most luminous ULXs (Lx > 10'^'^ erg/s) can be explained 
with a modest value of r/Edd = 0.1, favoring a larger IMBH 
mass. The enhancement of IMBH density around galactic 
disks by baryon compression may further increase the number 
of ULXs, possibly allowing a lower IMBH mass (see If 
the radiative feedback and ionization of ambient gas are sig- 
inificant, Cj should be raised to ~ 10 km s"' and the number 
of ULX would be significantly reduced, especially in the case 
of Mbh = 1OOM0. However, in the next section we will argue 
that this effect is not necessarily strong. 

6. ON THE RADIATIVE FEEDBACK 

The ULX luminosity is sufficient to ionize the surrounding 
medium; the Stromgren radius is 

rs - 7.0 X 10^%^'\[{^L%^ cm , (15) 

where e is the fraction of ionizing luminosity above 13.6 eV 
compared with Lx, and e_i = e/0.1. This is in most cases 
larger than the Bondi accretion radius, 

^ 9MEL = [J X lO^^M2vf cm. (16) 

Therefore, when an ULX ionizes the surrounding medium, 
the Bondi accretion rate should reduce, because of the in- 
creased C5 and expansion of gas driven by increased pres- 
sure. However, the gas that has already accreted to form an 
accretion disk will serve as the energy source of the ULX 
phenomenon for some time after the ionization. Based on 
the interstellar density gradient and angular momentum con- 
servation, the interstellar gas captured by an IMBH should 
form a rotationally supported disk at a radius ~ 7.5 x 

IQ^^mI^\~^°^^ cm (Agol & Kamionkowski 2002), which is 
much smaller than the Bondi accretion radius, r^- Assuming 
that the further accretion proceeds by a-viscousty in the stan- 
dard thin disk, accretion rate for ULX-level luminosity is not 
achieved until a sufficient amount of mass is accumulated and 
the surface mass density E at r ~ r^/ reaches that of the stan- 
dard disk for the mass accretion rate of ULXs, m. It would 



take a time of ^ InYir^ /m. Then ULX activity will last for 
a time scale of the disk mass consumption, InYir^/m, even 
after the ULX ionizes the surrounding medium. 

Therefore both the time scale of turning on of ULXs and 
duration of the ULX phenomenon are given by the accretion 
time scale of the standard disk. 



m 
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(The outer region of the three regions in the standard disk the- 
ory applies at r ^ rj.) It is interesting that this time scale for 
an ULX with Lx ^ lO'"' ergs s~' is similar to the crossing time 
scale fcross ^ lO*" yr of a n 10^ cm"-' cloud (^. This means 
that the duty cycle as an ULX could effectively be of order 
unity for an IMBH passing through a molecular cloud, even if 
the radiative feedback is significant. 

7. DISCUSSION AND CONCLUSIONS 

We conclude that the observed number of ULXs can be 
explained by halo population IMBHs passing through dense 
molecular clouds, which were produced as Pop III star rem- 
nants, if the Pop III star formation was very efficient as re- 
cently indicated by the unaccounted excess of the cosmic 
near-infrared background radiation. Luminosity function was 
calculated, and the ULXs of our scenario contribute to the ob- 
served luminosity function especially at the largest luminos- 
ity range, without overproducing sources in low-luminosity 
range where stellar-mass black holes are the dominant contri- 
bution. The ULX phenomena could last for 10^"'' yr and the 
duty cycle in a molecular cloud could be of order unity. 

We considered molecular clouds with a denstiy range of 
10^-10^ cm"^, but rather low-density clouds of « ~ 10^ cm"^ 
have the largest contribution to the expected ULX number It 
is consistent with observations of ULX locations; many ULXs 
are found in star forming regions but in outskirts of large con- 
centrations of gas (Zezas et al. 2002). 

Large ionized, expanding optical nebulae with a size of > 
100 pc are often associated with ULXs, with estimated total 
expansion energy of ^ 10^^ erg and age of ^ 10^ yr (PakuU 
& Mirioni 2002; Mushotzky 2004). Our estimate of the ULX 
duration is consistent with this age, and the total energy emit- 
ted by an ULX with Lx ^ lO'*" ergs s"' could be as large as 
~ 3 X 10^^ erg. A part of this energy may be emitted in 
UV and soft X-ray bands to ionize and expand the ambient 
medium, giving an explanation for the observed optical nebu- 
lae. 

The estimated ULX duration is not very different from the 
Eddington time scale, fedd = Mbh /{LEddMc^) = 4.4 x 10^?7-i 
yr. Therefore, some ULXs may have a chance to accrete gas 
for a time scale comparable to or longer than fEdd, by passing 
through a large molecular cloud or many molecular clouds. 
Such IMBHs will increase their mass considerably. In fact, 
the IMBH mass used in our calculation is not necessarily that 
of the Pop III remnants. The mass function of IMBHs in the 
local universe may have been altered from the original one at 
the formation epoch, by merging or gas accretion in the hi- 
erarchical structure formation. Such calculations have been 
performed by Islam, Taylor, & Silk (2004) and Volonteri & 
Perna (2005). According to these calculations, typical mass 
of free-floating IMBHs may be increased up to ~ IO'^Mq de- 
pending on model parameters, even if the initial mass of Pop 
III remnant BHs is less than 1000 Mq. The very large BH 
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mass suggested for the extremely bright ULX in M82 (Mat- 
sumo to et al. 2001; Kaaret et al. 2001) may be achieved by 
such processes. 

We used the popular NFW density profile for the Galactic 
halo, but the dark matter density could be significantly altered 
and increased around the region of the Galactic disk, by adi- 
abatic compression during the infall of baryonic gas to form 
stars. The number of ULXs would then be increased, because 
more IMBHs exist around the Galactic disk and the weighted 
mean surface density of molecular gas is also increased. Ac- 
cording to the model by Klypin et al. (2002), we estimate the 
boost factor of the ULX number as ^ 3. On the other hand, the 
same authors also shown that a modest exchange of angular 
momentum between baryonic and dark matter (through e.g., a 
rotating bar-like bulge) may result in almost no enhancement, 
or even decrease of the dark matter density. This should be 
taken into consideration as an uncertainty in the prediction 
made here, though the main conclusions are unlikely to be 
changed. 

Finally we give a list of the predictions for ULXs pro- 
duced by the scenario proposed here. Since ULXs are iso- 
lated IMBHs, optical counterparts (donor stars) will not be 



found. Future observations should confirm the trend that 
ULXs are preferentially found in star forming regions, but 
at outskirts of dense molecular clouds where the density is 
about 100 cm"^. Expanding nebulae with kinetic energy 
of ~ 10^^ erg should often be associated, without evidence 
for a past hypernova or multiple supernovae. The predicted 
slope of the IMBH luminosity function is a = 1 .5 at low lumi- 
nosity range, which is similar to that of HMXBs. Therefore, 
in a wide range of luminosity, we expect that a small frac- 
tion of low-luminosity X-ray sources are the same population 
with IMBH-ULXs, showing RJAF spectra without compan- 
ion stars. It is interesting to search for such a population in 
the Galaxy or nearby galaxies. The excess of CNIB must be 
confirmed with greater significance. The high density of free- 
floating IMBHs in galactic haloes may be detected by future 
gravitational lensing experiments (e.g., Totani 2003). 
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Fig. 1. — Left: The X-ray luminosity function of free-floating IMBHs (Mbh =100 Mq) passing through dense molecular clouds, normalized to SFR>5 = 
5OM0 yr~'. For each line marking, there are two model curves: the upper one for rjEdd = 1 and the lower one for = 0.1. The sound speed is set as Cs = 0.3 
and 10 km s~' for each thick and thin curve, respectively. The dynamical friction parameter is set to be ^ = 1 and 10 for each soUd and dashed curve, respectively. 
The data points show the observed universal X-ray luminosity function from Grimm et al. (2003). Right: Same as the left panel, but for Mbh = IO'M© . 



